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Abstract 

Two widely proposed fcf-dependent gluon distributions in the small- x saturation regime are 
investigated using two particle back-to-back correlations in high energy scattering processes. The 
Weizsacker- Williams gluon distribution, interpreted as the number density of gluon inside nucleus, 
is studied in the quark-antiquark jet correlation in deep inelastic scattering. On the other hand, the 
unintegrated gluon distribution, defined as the Fourier transform of the color-dipole cross section, 
is probed in the direct photon-jet correlation in pA collisions. Dijet-cor relation in pA collisions 
depends on both gluon distributions through combination and convolution in the large N c limit. 
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As the foundation of high energy hadronic physics, QCD factorization enables us to sepa- 
rate the short distance perturbative physics from the long distance non-perturbative effects. 
Its prediction power relies on the universality of the parton distributions among different 
processes. Recent studies [11-01 have shown that the naive fc t -factorization is violated in two- 
particle production in hadron-hadron collisions. Nevertheless, in this letter, we establish 
an effective factorization in hard processes in nuclei scattered by a dilute probe by modi- 
fying the parton distributions at small-x. Although the parton distributions are no longer 
universal, they can be constructed from several universal individual building blocks. 

The saturation phenomena in high energy hadronic interactions have attracted great 
attention in recent years. It has long been reco gniz ed that the gluon dynamics in QCD 



at small- x is responsible for these phenomena |8Hlll]. An effective theory the color glass 
condensate (CGC), was proposed to systematically study this physics On the exper- 
imental side, there exist strong indications of the saturation from the structure function 
measurements in deep inelastic scattering (DIS) at HERA and the shadowing effects found 
in inclusive hadron production in dA collisions at RHIC [llj. Ongoing and future experi- 
ments from both RHIC and LHC, and in particular, the planned Electron- Ion Collider [121 ]. 
shall provide further information on this. 

An important aspect of the small-x gluon distribution function in nucleons and nuclei 
is the resummation of the multiple interactions of the hadronic probe with the target, be- 
cause the gluon density is so high and these interactions have to be taken into account. In 
order to study the resummation effects, a transverse momentum dependence is introduced. 
They are referred to as the fc t -dependent gluon distributions, also called unintegrated gluon 
distribution (UGD) functions. Two different forms of the UGDs have been used in the lit- 
erature. The first gluon distribution, also known as the Weizsacker- Williams (WW) gluon 



distribution, measures the number density of gluons in the CGC formalism 10( , whereas the 



second one defined as the Fourier transform of the color dipole cross section, appears in the 



calculations for, e.g., single inclusive particle production in pA collisions [11]. However, it 
has been argued that they can not be distinguished especially in the small k± region though 
they differ dramatically, and both of them are often used jl3| . 

In this paper, we study two particle correlations in various high energy scattering pro- 
cesses as probes to these UGDs. There have been intensive investigations of these processes 
in the last few years [l|-0]. Taking the quark distribution as an example, it was shown in 
Ref. [7[ that, in the large nuclear number limit, an effective factorization can be achieved 
with modified parton distributions of nucleus in pA and 7*A scattering processes where 
the multiple interaction effects can be resummed in the small-x formalism. Following this 
argument, we focus on the processes with a dilute system scattering on a dense target, 

B + A^H x {k x )+H 2 (k 2 )+X , (1) 

where A represents the dense target (such as a nucleus), B stands for the dilute projectile 
(such as nucleon or photon), H\ and H 2 are the two final state particles with momenta 
ki and k 2 , respectively. We are interested in the kinematic region where the transverse 
momentum imbalance between them is much smaller than the individual momentum: q± = 
\ki± + k 2 ±\ <C P± where P± is defined as (ki± — k 2 ±)/2. This is referred as the back-to- 
back correlation limit (the correlation limit) in the following discussions. An important 
advantage of taking this limit is that we can apply the power counting method to obtain the 
leading order contribution of q±/P± where the differential cross section directly depends on 
the UGDs of the nuclei. For example, the quark- ant iquark jet correlation in deep inelastic 
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scattering (DIS) directly probes the first type of the UGD, whereas the direct photon- 
quark jet correlation in pA collisions probes the second type of UGD. The dijet (di-hadron) 
correlations in pA collisions can probe both gluon distributions, though the connection is 
more complicated. 

Let us first discuss the conventional gluon distribution, generalized to include the trans- 



verse momentum dependence [14j, [15 

xG {1) {x,k ± ) -- 



-e 



x(P\F + *(C^±)CtCoF +t (0)\P) , (2) 



(2vr) 3 P 

where F^ v is the gauge field strength tensor F£ u = d^A v a — d u A^—gf abc A^A v c with f abc the an- 
tisymmetric structure constants for SU (3), and C^ — V exp{—ig J^°? d(~A + ((, £±)}V exp{ — ig d(j 
Aj_((~ = oo, is the gauge link in the adjoint representation A^ = A^t a with t a = —ifabc- 
It contains a transverse gauge link at spatial infinity which is important to make the def- 
inition gauge invariant |16| . These gauge links have to be made non-light-like to regulate 
the light-cone singularities when gluon radiation contributions are taken into account (lij . 
This does not affect the following analysis, because there is no light-cone singularity in 
the calculation. By choosing the light-cone gauge with certain boundary condition for the 
gauge potential (for example, in the above definition, A±((~ = oo) = 0), we can drop out 
the gauge link contribution, and find that this gluon distribution has the number density 
interpretation. Then, it can be calculated from the wave functions or the WW field of the 



nucleus target [10|, [17|. At small- a: for a large nucleus, it was found 10 



wiw i \ S± N 2 -l f d 2 r± e~ ik± ' r± i 

xGW{x,k ± )= / c / — — 1-e-— , 3 

7T 2 a s N c J (2TT) 2 rj_ 



where N c = 3 is the number of colors and Q s is the gluon saturation scale [llj. We have 



cross checked this result by directly calculating the gluon distribution function in Eq. (E 



following the similar calculation for the quark in Ref. [16|, [18 . 

Despite the nice physical interpretation, it has been argued that the gluon distribution 
in Eq. fl5]) is not directly related to physical observables in the CGC formalism. However, 
we will show that xG^ can be directly probed through the quark-antiquark jet correlation 
in DIS, 

j* T A q(h) + q(k 2 ) + X . (4) 

The calculations are performed for Q 2 in the same order of Pj_ and the quarks are massless. 
Extension to the real photon scattering and/or massive quarks in the final state is straight- 
forward. We will show the results for the transversely polarized photon, and that for the 
longitudinal one follows accordingly. We plot the typical Feynman diagram for the process 
of (jlj) in Fig. 1, where the bubble in the partonic part represents the hard interaction vertex 
including gluon attachments to both quark and antiquark lines. Fig. 1(a) is the leading 
Born diagram whose contributions can be associated with the hard partonic cross section 
times the gluon distribution from Eq. ([2]) Q. In high energy scattering with the nucleus 
target, additional gluon attachments are important and we have to resum these contribu- 
tions in the large nuclear number limit. Figs. l(b,c) represent the diagrams contributing at 
two-gluon exchange order, where the second gluon can attach to either the quark line or the 
antiquark line. By applying the power counting method in the correlation limit (q± P±), 
we can simplify the scattering amplitudes with the Eikonal approximation j^J. For example, 
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(d) (e) (/) 

FIG. 1. Typical Feynman diagrams contributing to the quark-antiquark jet correlation in deep 
inelastic scattering (a,b,c) and photon-jet correlation in pA collisions (d,e,f): (a,d) leading order, 
where the bubble represents the gluon attachments to both quark lines; (b,c,e,f ) two-gluon exchange 
diagrams. 



Fig. 1(b) can be reduced to: gj^—q^ + ie)T b T a in the above limit, where q 2 is the gluon mo- 
mentum, T b is the SU (3) color matrix in the fundamental representation and r a represents 
the rest of the partonic scattering amplitude with color indices for the two gluons a and b. 
Similarly, Fig. 1(c) can be reduced to: —g/{—qt + 2e)r a T b . The sum of these two diagrams 
will be g/{-q£ + it) [T b T a - T a T b ] . Because of the unique color index in T a , we find the 
effective vertex as, 

Fig. l(b,c) ±^(-i g )(-if bca )T c , (5) 

-qj + it 

which is exactly the first order expansion of the gauge link contribution in the gluon distri- 
bution defined in Eq. (j2J). Here, the crossing diagrams of Figs.l(b,c), by interchanging the 
two gluons, are included in the calculations to compare to the gauge link expansion results 
(see also Ref. pl|). For all high order contributions, we can follow the procedure outlined 



in Ref. [16] to derive the gluon distribution. Therefore, we obtain the following differential 



cross section for the quark-antiquark jet correlation in DIS process 

da7 *A^qq+X 



dV.S. 



S(x r - l)x g G (1) (xg, q±)H 7i , g ^qq , (6) 



where x g is the momentum fraction carried by the gluon and is determined by the kinematics, 
£ 7 » = z q + Zq with z q and z q being the momentum fractions of the virtual photon carried 
by the quark and antiquark, respectively. The phase space factor is defined as dV.S. = 
dyidy 2 d 2 P±d 2 q±, and y\ and y 2 are rapidities of the two outgoing particles. The leading 
order hard partonic cross section reads H 7 ^ g ^. qq - = a s a em e 2 q (s 2 + Q A ) j '(s + Q 2 ) A x [u/t + t/u) 
with the usually defined partonic Mandelstam variables s, t and u. By taking Q 2 = 0, we 
can extend the above result to the case of dijet productions in real photons scattering on 
nuclei. The above process (j4]) can also be analyzed following the procedure in Ref. 0, where 
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the gluon distribution Eq. (j2J) is written in the fundamental representation, 



xG (1 \x, k±) 



_ e ixP+£, -ifc ± -£i 



x (P|Tr [F +i (r, £±)U [+ ^F +l (0)U [+] ] |P) , (7) 

where the gauge link = U n [0, +00; 0] U n [+00, £J with f/ n being the light-like Wil- 
son line in covariant gauge. 

The quark- ant iquark jet correlation in DIS (JJJ) can also be calculated directly in the color- 



dipole and CGC formalism [201 ] . The basic element is the qq (color-dipole) wave function of 
the virtual photon, combined with the multi-scattering of the dipole on the nuclear target. 
Following this formalism, the amplitude can be written as 

\ A \2 _ N 2 f d 2 £ dV d 2 b d 2 b' _ ifelx . (a ._ a ., ) 
lAl - 7VcCW3 9 J (27r) 2 (27r) 2 (27r) 2 (27r) 2e 
Xe -ik 2X -(b-b>) J2i>r( x - b)ip T {x' - b') 

x \ + S${z, b; b>, x') - S<g(x, b) - S®{1/, x 1 )] , (8) 

where ip? is the qq Fock state wave function of a transversely polarized photon depend- 

ing on e 2 = z(l - z)Q 2 with z = z q , Sjg{x,b) = ^ (TrU{x) W{b) ) v 5^ (x, 6; 6', x') = 

(TrU (x)W \x')V (b')U^(b)) x , and U is the Wilson line describing the multiple scatterings 
of a single quark with the nuclear target. The expectation value of multiple Wilson lines can 



be found in Refs. |2CH22| . In order to study the differential cross section in the correlation 
limit, we substitute u = x — b and v = zx + (1 — z)b. Thus, the exponential factor becomes 
e -iqx-(v-v ) e -tPx-(«-M ) w here P± = (1 — z)k\± — zk2± ~ P±. Then, we can expand the in- 
teraction part of the bracket at small u and u', but keep v and v' fixed. We find that the 
remaining contribution comes from the term involving inelastic scatterings while the elastic 
scattering part cancels out. Therefore, the square bracket in the above equation becomes 
jfUiii'j (TrdiU(v)W (v')djU(v')W(v)^ x . With this expansion result, we further find that 

the wave function integral with UiUj leads to 5ij(Pj_ + £f)/{P± + e/) 2 , an d the differential 
cross section can be simplified as, 



da -y^A^qqX p4 + A 

dvs = a em e 2 q a s 5 (x r - 1) z(l - z) (z 2 + (1 - z) 2 ) ^ + e 2y 

X (l 67r 3 ) f <^^ e -n^v-v') 2( ^ F i+r v)umF i + (g) 

J (27r) b x 9 

To compare with the differential cross section in Eq. ([6]), we notice that the partonic 
Mandelstam variables can be expressed in terms of Pj_ and z as: s — P]_/(z(l — z)), 
i = —(P 2 +e 2 )/(l — z), and u = —(P]_ + e 2 )/z. Substituting these relations into Eq. (9), tak- 
ing the small-x approximation for the gluon distribution, and correcting the normalization 
for the states in the calculation of the associated matrix elements, we find that it agrees with 
the factorization result ([6]) completely. This consistency is very important to demonstrate 
that we have a unified picture for the quark-antiquark correlation in DIS. 
On the other hand, the direct photon-quark jet correlation in pA collisions, 

P A -> 7(h) + q(k 2 ) + X , (10) 
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probes a different gluon distribution. We plot the relevant diagrams in Fig. 1 (d,e,f), again 
for the leading one gluon exchange and two gluon exchanges. Similarly, the two gluon 
exchange contributions can be summarized as 

Fig. l(e, f) ^rr-i-ig) {T b T a + T a T b ) , (11) 

where the plus sign comes from the fact that the second gluon attaches to the quark line in 
the initial and final states. Since there is no color structure corresponding to Eq. (fTTj) . we 
can only express it in the fundamental representation. Following Ref. 0, we find that the 
gluon distribution in this process can be written as 

xG (2) (x,k ± ) = 2 J ^fp + elXP+r ~ lk±<± ( P \ Tl [F +l (C^±)U H] F +t (0)U [+] ] \P) , (12) 

where the gauge link = U n [0, — oo; 0] U n [— oo, This gluon distribution can also 

be calculated in the CGC formalism where it is found to be xG^(x,q±) ~ S±F Xg (q±) 
with the normalized unintegrated gluon distribution F Xg (q±) = J ^pe~ iq - L ' rj -S x 2 g > (0, r±). 
Therefore, the differential cross section of (TT0|) can be written as 

— ^ = 2_^x 1 q(x 1 )x g G {2 \x g ,q ± )H qg ^ q , (13) 

/ 

where q(xi) is the integrated quark distribution from the projectile nucleon. Because we are 
taking large nuclear number limit, the intrinsic transverse momentum associated with it can 
be neglected compared to that from the gluon distribution of nucleus. The hard partonic 
cross section is H qg ^. iq = a. q a P e 2 / (N r s 2 ) (—s/u — u/s). We can calculate process (flOl) in the 
CGC formalism directly [23j], and again we find that these two calculations are consistent 
with each other. 

The gluon distributions defined in Eqs. ([2]) and f|T2|) have the same perturbative behavior 
and they reduce to the same gluon distribution after integrating over k±. However, they do 



differ in the low transverse momentum region [ll|, |13| . It will be very important to measure 



the low k± behavior of the correlation processes of (j3J) and ffTUl) to test these predictions. In 
particular, the planed EIC machine will be able to study the quark-antiquark jet correlation 
in DIS process, whereas RHIC and future LHC experiments shall provide information on 
direct photon-quark jet correlation in pA collisions. 

The extension to the dijet-correlation in pA collisions is straightforward. The relevant 
initial and final state interaction phases have been calculated in Ref. 0]. To obtain the 
differential cross section in the correlation limit, we have to take the large iV c limit and the 
mean field approximation 0, 24j. After a lengthy calculation, we find 



^ a (pA^Di}et+X) 

dV.S. 

= y]Xiq(xi)-^ [-Fqg HqgUqg + Fqg 1 Hqg-tqg] 



a 2 



a 2 

+ x 1 g(x 1 )-^ 



jr(i) ( _ _i_ 



+Jr (2) / R m _ , l H (2) \ , lr(3) ff (3) 
J gg \ 11 gg^qg ~ 2 99^99 J ^2" 99^-gg 



(14) 
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where again q{x\) and g(xi) are integrated quark and gluon distributions from the projectile 
nucleon. The hard partonic cross sections are defined as 



(1) u 2 (s 2 + u 2 ) (2) _ s 2 (s 2 + u 2 ) 
*w -2§ut 2 ' q9 ^ qg -2sut 2 
H (D 1 2(t 2 + £ 2 ) 2 r{2) _ 1 4(t 2 + u 2 ) 



>gg-*n ANc §2{li > "ss-> M 4jVc S 2 
F (i) _ 2p + fl 2 )(g 2 -tfl) 2 ^ (2) _4(S 2 -t«) 



~<3) 2 ( g2 -^)" 



M 2 t 2 S 2 ' 9 ™ UtS 2 

, 2 

(15) 

and the various gluon distributions of nucleus A are defined as 

J-W = *G( 2 >(a;,g ± ), jg> = J xG^ ® F (g 2 ) , 

= / zG (2) (gi) ® F( fe ), Jg) = - / ^^xG( 2 )( gi ) ® F(g 2 ) , 

J*f = j xG^( qi ) ® F(g 2 ) ® F(g 3 ) , (16) 

where (g) represents the convolution in momentum space: f <g) = J d 2 q 1 d 2 q25^(q± — qi — 
q-i). Clearly, this process depends on both UGDs in a complicated way, and the naive k t - 
factorization does not hold. We have checked the above results in different kinematics. First, 
we recover the inclusive dijet cross section by integrating over q±. Second, in the dilute limit 
of A, or equivalently at large q±: P± 3> _q± 3> Qs,^-qcd, we reproduce the dijet-correlation 
in the collinear factorization approach |3(. Last but not least, we find the CGC calculation 
agrees with above results perfectly. 

Recently, both STAR and PHENIX Collaborations have published experimental results 
on di-hadron correlations in dAu collisions, where a strong back-to-back de-correlation of 
the two hadrons was found in the forward rapidity region of the deuteron [25| . These 
results have been compared to a number of theoretical calculations in the CGC formalism, 



where different assumptions have been made in the formulations [26|-|28|. In particular, the 
numerical evaluation in Ref. [26[ only contains the first term in the qg channel in Eq.(14). 
The missing term is equally important and should be taken into account to interpret the 
STAR data. We plan to compare our results to these calculations and those in Ref. [29| for 
qq production in pA collisions in a future publication [13], together with detailed derivations 
of generalized results illustrated in this paper. 

In summary, we have studied the fc t -factorization for dijet production at small-x in nuclei. 
We found that different gluon distributions probed in different dijet production processes 
can be built from two basic building blocks, the Weizsacker- Williams distribution and the 
Fourier transform of the dipole cross section. The most important result is that the DIS 
dijet process can directly measure the well-known WW gluon distribution which would be 
very interesting physics to pursue at the planned EIC. 
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